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ABSTRACT 


A  priori  estimates  for  solutions  of  the  quasilinear  hyperbolic-parabolic 
equations  governing  the  initial  value  problem  describing  the  motion  of 
compressible,  viscous  and  heat-conductive,  Newtonian  fluids  are  derived  by 
means  of  a  new  energy  method.  This  technique  enables  us  to  simplify  and  unify 
our  previous  results  on  the  global  existence  in  time  and  uniqueness  of  smooth 
solutions  of  these  equations  for  sufficiently  smooth  and  "small"  initial  data 
and  to  obtain  their  rate  of  decay. 
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SIGNiriCANCE  AND  EXPLANATION 


The  motion  of  compressible,  viscous  and  heat-conductive  Newtonian  fluids 
is  described  by  a  system  of  partial  differential  equations  which  is  of 
hyperbolic-parabolic  type  and  highly  nonlinear.  One  of  the  first  mathematical 
problems  associated  with  this  system  is  the  initial  value  problem.  We  obtain 
the  existence  of  a  unique  smooth  global  solution  in  time  for  the  initial  value 
problem  and  also  the  decay  rate  of  the  solution  as  time  tends  to  infinity. 
Since  the  system  is  quasilinear  with  respect  to  the  unknowns:  density, 
velocity  and  temperature,  we  need  to  assume  that  the  initial  data  are  close  to 
the  constant  equilibrium  state.  The  purpose  of  this  paper  is  to  obtain  a 
priori  estimates  for  the  solutions  of  these  equations  by  means  of  a  new  energy 
method.  This  technique,  although  still  necessarily  laborious,  enables  us  to 
simplify  and  unify  our  previous  results,  described  briefly  in  the  abstract  and 
obtained  jointly  with  T.  Nishida  (see,  e.g.  MRC  TSR  »1h91). 
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The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


AN  ENERGY  METHOD  FOR  THE  EQUATIONS  OF  MOTION 
OF  COMPRESSIBLE  VISCOUS  AND  HEAT-CONDUCTIVE  FLUIDS 


Akitaka  Matsumura 

^ 1 .  Introduction  and  Main  Theorem. 

In  previous  papers  [1],  [2],  we  have  investigated  the  global  solution  in 
time  of  the  initial  value  problem  for  the  following  equations  governing  the 
motion  of  isotropic  Newtonian  fluids; 


(1.1) 


+  (pu  -  0  , 

j 


u^  +  u^u^  +  ■“  P  =  —(yCu^  +  u^  )  +  p'u*^  i  “  1,2,3, 

t  X,  pPx,  p^»^  X,  X,  x^  V 


0p 


+  u^9  +  — ^  u^  =—  (k6  )  +4')  , 


’‘j  <^V  ’‘j  P^V 


X  .  X  . 

j  D 


with  the  initial  data 


(1.2) 


(p,u,0)(0,x)  =  (  Pp  ,Ujj,  )  (x)  , 


3  12  3 

where  t  >  0,  x  =  (x^,X2,x^)  e  R  ,  p  is  the  density,  u  =  (u  ,u  ,u  )  is  the 
velocity,  9  is  the  absolute  temperature,  p  =  p(p,0)  is  the  pressure, 
u  =  ii(p,9)  and  p'(p,9)  are  viscous  coefficients,  <  =  ic(p,9)  is  the 
coefficient  of  heat  conduction,  c^  =  c^(p,9)  is  the  heat  capacity  at 
constant  volume,  and  4'  =  ■-  (u^  +  u'*‘  )^  +  u'  (u^  is  dissipation 

^  X,  X  .  X  , 

K  1  1 
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function.  We  shall  consider  the  solutions  only  In  a  nelqhbourhood  of  any 
fixed  constant  state  (p,u,9)  »  (p,0,0)  where  p,0  are  any  positive 
constants.  Moreover,  we  shall  make  the  following  natural  assumptions  on  the 
hyperbolic-parabolic  system  (1.1)  throughout  this  paper; 

(1)  p,  c^,  u,  u'  and  <  are  smooth  functions  of  0  =•  {(p,u,0)  ;  (p  -  p|, 

|u|  ,  I  0  -  0|  <  }  , 

2 

(11)  p  ,  p  ,  c  ,  u»  <  >  0  and  p*  +  —  u  >  0  in  0, 
p  0  V  3 

where  <  min(p,9). 

In  [1],  we  succeeded  in  obtaining  a  global  solution  in  time  of  (1.1), 
(1.2)  by  using  energy  methods  which  were  rather  technical  and  complicated 
under  the  assumptions  that  the  fluid  is  an  ideal  and  polytropic  gas  and  that 

(Pq  -  p»  Up,  -  0)  is  suitably  small  in  , 

V 

where  H  represents  the  usual  Sobolev's  spaces  with  the  norm  also 

proved  the  decay  of  the  solution  to  the  constant  state  ( p,  0,  0),  but  we 
were  not  able  to  estimate  the  decay  rate  of  the  solution.  In  [2] ,  we  obtained 
both  the  global  solution  in  time  of  the  original  problem  (1.1),  (1.2)  and  its 
decay  rate  to  the  constant  state  by  using  a  calculation  of  the  decay  rate  for 
the  linearized  equations,  together  with  the  energy  estimates.  Moreover,  we 
had  to  investigate  the  precise  properties  of  the  spectrum  of  linearized 
equations,  and  assume  that 

(pp  -  p,  Up,  Bp  -  9)  is  suitably  small  in  n 


The  purpose  of  this  paper  is  to  employ  a  different  energy  method  to 
handle  the  nonlinearity.  The  present  approach  is  simpler  and  less  technical 
than  the  energy  form  used  in  [11,  moreover  we  do  not  require  all  of  the 
precise  properties  of  the  linearized  equations  obtained  laboriously  in  [2] . 
However,  we  note  that  because  of  the  relative  roughness  of  the  new  method,  the 
coefficients  of  the  various  estimates  might  be  more  rough  than  those  in  [1] 
and  [2].  By  making  use  of  this  approach,  we  obtain  both  the  existence  of 
global  smooth  solution  in  time  as  well  as  its  decay  rate  to  the  constant  state 
for  the  general  case  (1.1),  (1.2)  under  the  assumption  that 

—  —  3 

(p^  -  p,  u^,  -  9)  is  suitably  small  in  H  , 

(the  previous  approach  required  n  L^).  This  method  will  be  applied  to  an 
initial  boundary  value  problem  for  (1.1),  (1.2)  in  a  subsequent  paper. 

The  main  result  is; 

Theorem  1.1.  Consider  the  initial  value  problem  (1.1),  (1.2)  and  suppose 
the  initial  data  ( p^  -  p,  u^,  -  8)  e  H^.  Then  there  exist  positive 

constants  and  Cg  such  that  if  (p^^  -  p,  u^,  0^  -  911^  <  Eg,  the  problem 

(1.1),  (1.2)  has  the  unique  global  solution  in  time  satisfying 

p  -  p  e  C°(0,  +«>;  H^)  n  c\o,  +c»;  H^)  , 

(u,  9  -  9)  e  c'’(0,  +<»;  H^)  n  c’(0,  +'»;  h’)  , 

and 

3 

sup  |(p-p,  u,  9  -  0)(t)(  <  C^(1+t)  ^  ilPi^-p,  ,  0^  -  011^ 


(1.3) 


X 


52.  Proof  of  Theorem 


First  rewrite  the  problem  (1.1),  (1.2)  by  the  chanqes  of  variables: 
p+p+p,  u+u,  0+0+9,  p( p+p,  6+0)  +  p(p,9),  m(p+p,  0+6)  +  y(p,9) 

so  on  resulting  in 


p^  +  (p  +  p)u^  +  u^p^  -  0  , 

j  j 


i  i  i  ”  6 

u^+u^u  +PP  +  -  9 

t  X,  p  X.  —  X, 

]  1  p+p  i 


(2 


... 


—  (u(u^  +  )  +  u'u’^  6^^) 

—  X  .  X.  X,  X  , 

p+p  :j  i  X  -J 


1  4  1 

0  +  +  — - 2  yj  - -  )  +  f)  , 

’‘j  (P+P)c  ’‘j  (p+p)c  ’'j  ’‘j 


(2.2) 


(p,u,9)(0,x)  =•  (Pq,  Uq,  SqXx) 


Define  a  positive  constant  bv  the  Sobolev's  lemma  so  that  for 

ilfll^  <  Eg  we  have  supifl  <  Cllfll^  <  Yg*  Denote 
.  QL  a. 

n'^f  *  {3  f/3x^  Sx^  3x^  for  all  a,  +  Oj  +  =  k)  and  define 

11*11  »  "^ben  the  solution  of  (1.1),  (1.2)  is  sought  in  the  set  of 

functions  X(0,  +<*’;  E)  for  some  E  <  E^^,  where  for  0  <  t^  ^  '"2  ^ 
x(t,,  tjj  E)  -  {(p,u,9)  :  p  e  C°(t^,  t^;  H^)  n  c\t^,  t^;  H^)  , 

n„  c  L^(t^,  t^;  H^),  (1  p  e  L^(t^,  t^;  L^)  , 

( 1+t)D^p  G  L^(t^,  t,;  L^)  , 

(u,1)  c  t^;  H^)  n  c\t^,  t^;  h’)  , 

r)(u,0)  G  L“{t^,  t^;  H^),  ( i+t)'^2Q2j^^gj  e  L^(t^,  t^i  l^)  , 

(  1+t)n\u,  0)  G  L^(t^  ,  t^;  h’  )  , 

and  N(t,,  t^)  <  F,  (E  <  H  • 


and 


We  define  N(t^,  tj)  by 

N^(t  ,  t  )  »  sup  ( II  ( p»u,  6)  ( t)  llj  +  t  IID(  p,u,  0)  (t)  • 

ti<t<t2 

.  IID^(p,u,0)(t)  B^)  +  IIDp(T)ll2  +  I1D(u,0)(t)  + 

+  T  110^  (  p/U,  0)  (  t)  +  T^(  IID^p(  t)  +  BD^  (u,  0)  (  t)  B^)  dT 

We  prove  Theorem  1.1  by  a  combination  of  a  local  existence  result  and  a  priori 
estimate  for  the  solution  in  X, 

Theorem  2.1.  (local  existence)  Consider  the  initial  value  problem  (2.1) 
for  t  >  t^  with  the  initial  data  at  t  »  t^  ^ 

(2.3)  (p,u,0)(t^)  e  . 

Then  there  exist  positive  constants  e^,  C^(e^C^  <  E^^)  and  t  which  are 
independent  of  t^  such  that  if  N(t^*  t^)  <  »  the  problem  (2.1)-(2.3)  has 

the  unicTue  solution 

(p,u,0)  e  X(t^,  t^  +  Tf  C^N(t^,t^))  . 

The  proof  of  Theorem  2.1  is  given  in  the  same  way  as  in  [1]. 

Theorem  2.2.  (a  priori  estimates)  Suppose  that  the  initial  value 
problem  (2.1),  (2.2)  has  a  solution 

(p,u,0)  e  X(0,T;  E) 


-5- 


for  some  T  >  0  and  some  E  <  E^.  Then  there  exist  three  positive  constants 

/  ^ 

^2  ^3  ^  ^  ^0^  which  are  independent  of 

T  such  that  if  "Po''^0'^o"3  ^  ^2  E  <  e^,  then  the  solution  satisfies 

the  a  priori  estimate 

(p,u,9)  e  X(0,T;  11^)  . 

Proof  of  Theorem  2.2.  Take 


•  r  S  ^1 


We  may  use  the  standard  continuation  arguments  of  a  local  solution  on 

[0,nT],  n  =  1,2,...  to  get  the  global  solution.  In  fact,  by  the  local 

existence  theorem,  the  definition  of  and  the  assumption 

lip  ,,u„,9„ll,  <  E^,  we  have  a  local  solution 
■003  0 


(p,u,0)  e  X(0,T7  '■^O  '  3  ^ 


Ry  ilp^,  C  and  llpp  ,U|^  ,  llj  <  C^'  ®  priori  estimate 

gives 


(p,u,6)  E  X(0,t;  ^2"^0''^0 


®0"3’ 


Then  bv  ^2 '' ^0  ' '^0  '  ^0  Vl  ^  ^2^0  ^  local  existence  theorem  with 

T,  we  have  again 


(  0,u,  0)  E  X( 
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Noting  that 


N^(0,2t)  <  N^{0,t)  +  N^(t,2t)  , 


we  also  have 


(p,u,0)  e  X(0,2t»  Cj/I+C^  "Po'’^0'®o"3'  * 


/  2  7  2 

Now  by  Cj/I+C^  *^0''^0'®0*3  ^  ^0  ^  ^3'  *  priori  estimate  shows 


{p#u#0)  c  X(0«2ti  It  pQ !  ^^0^3^  * 


Thus  we  can  continue  to  use  the  same  arguments  on  [nT>  (n-t-Dil  and 
[0,  (n+1)T]  successively  n  «  2,3,...  .  Finally  the  estimate  (1.3)  follows 
from  Nlrenberg's  inequality  [3] 


(2.4) 


sup|p,u,0|  <  Clip, u, 011^  IID^  (  p,u,  0)  B** 


§3.  A  Priori  Estimates 

We  present  here  an  energy  method  to  obtain  a  priori  estimates  for  small 
solutions  of  equations  with  disslpatin.  First  we  rewrite  the  system  (2.1)  so 
that  all  the  nonlinear  terms  appear  at  the  right  hand  side  of  equations; 


-7- 


(3.1) 


L°  =  p  +  pu^  «  , 

t 

\  Pl^x.  ^  P2\  ■  ^x.x^ 


4  i  »4 

L  5  0  +  P,u  -  K0  =  t  , 

t  3  X .  x .X . 

3  3  3 


(  U+y’  )u^  =  f  ^  , 

XX 

i  3 


where 


=  p^(0,0)/p,  =  Pg{0,0)/p,  Pg  =  9pg(0,0)/p  C^(0,0)  , 


p  2 

U  =  y(0,0)/p,  U’  =  y*(0,0)/p,  <  =  <(0,0)/p  C  (0,0)  , 


=  -  pu^  -  u^p  , 

X.  X. 

3  3 

+  (p,  -  p  /(p+p)p  +  (P-  "  p  /(p+p))0  + 

X.  1  p  X,  2  a  X 

3  i  ^ 

+  (u  p  +  u„0  ){u^  +  )/(p+p)  + 

3  3  3  i 

+  (y'p  +  u'0  )u^  /(p+p)  +  (y/(p+p)  -  U)u^  + 

p  Xi  9  X.  X.  x.x. 


+  ( (u+u* )/(p+p)  -  (u+u’))u 


x.x. 
1  3 


an'l 

=  -  u’e  ^  (P3  -  (9^0)Pg/(p+p>c^)u^  +  (<^p^  + 

j  j  3  3 

+  (K/(p+p)C^  ■  ^’®X  X  ' 

3  j 

Define,  for  k,  Z,  m  =  0, 1,2,3, 


A’'(t)  =  /  ^  n’^f”  .  d’'p  +  D’'f^  .  +  —  n'^f'^  .  d'^o  dx  , 

P  ^3 


m  /  ni  0  vi\  o  ...W  J 

B  (t)  H  f  D  f  ‘Dp  +  ■  J"  D  f  ‘Dp  dx  , 

’'i  ’'i  2y+u’  ’'i 


•ft)  ;  /  — ^  o^'u^  •  dx  . 

2^+u’  ’‘1 


/ ( p+p)C^ 
3 
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Then  we  have  the  following: 

Lemma  3.1.  There  exist  positive  constants  v  and  C 
of  t  such  that 

(1)^  IID’^(p,u,e)  (t)  11^  /q  f 

<  C(  IId''(Pq,Ujj,0q)II^  +  Jj  lA’"(T)ldT), 

(ii)^  IID"‘p(t)ll^  /q  ^ 

<  CCKD^PqB^  +  «d"'  +  “D™  ^u(t)ll^ 

+  iid"'(u , 0) ( t)  11^  +  Ib'"  \t)I  +  Ic"* 

o  o  ^t.  P.  k'f  1  2 

(iii)„  ,  t  HD  (p,u,0)(t)ll  +  V  f  T  HD  (u,0)(t)H 

Ilk  0 

<  c( T*’  ^  1ID*^(p,u,0)(t)  +  I  t)  1  dx) 

(iv)  t^llD^^pi  t)  t*'1d"*p(  x)  ll^dx  < 

Z,m  "0 

<  C(t^llD'"'\(t)  11^  +  /p  x^”^  11d"'p(  x)  +  x^  ^ 
+  x^Hd'^(u,  0)  (  x)  +  1  X^c"  ^(x)l  +  1  x^B*" 

i  =  1,2,  k  =  2,3  . 


which  are  Independent 

0  <  k  <  3  , 

+ 

■\x)|dx),  1  <  m  <  3 

dx 

,  Z  =  1,2,  k  =  1,2,3 

m-1  ,  2 

HD  u(x)U  + 

\  x) I dx  , 
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Lemma  3.2.  There  exist  a  positive  constant  C  which  is  independent  of  t 
such  that 


(3.2)  N  (0,t)  <  Cdip  ,UQ,9gll2  +  /q  I  (1+t^)(|A^(t)|  +  |b^(t)|  + 

4=0 


+•  |C*'(t)I  )  +  (1  +  T^)|A^(T)|dT) 


I>roof  of  Lemma  3.1.  For  (i),,  and  (iii).  .  ,  we  may  estimate  the  equality 

^  Jt  /  K 


f  —  n''(L°)  •  d’^p  +  d''(L^)  •  +  —  d’'(L^)  •  D^'e  dxdT 

Jo  -  P3 


ft  Ik 

T  A  ( T)dT  , 


which  impli^^s  after  integratinq  by  parts 


i  ,  ^1  ,  k  ,2  1  ,  k  ,2  ^2  ,  k„,2, 

t  /—  iDpI  +-|Du|  +  - —  I D  0 1  dx 

2p  *^3  t=0 


rt  2,  f  k+1  ,2  .1  k  j  .2  '^^2  2.  k+1  , 

+  J-T  (uln  u|  +(g+ij')|Du'^|  +  -  T  |D  6 1 

'  P3 

ft  r  ,  k  ,2  1  ,  k  ,2  ^2  ,  k„,2, 

=  In  J,T  (1-6  )  J  —  Id  pi  +  -  Id  u|  +  - —  Id  9|  dx  + 

"  2;  2  2^3 


dxdi 


0  k 

+  T  A  ( T )dT  , 


where  6^'^  represents  Kronecker's  delta.  These  inequalities  prove  (i))^ 


and  (  !.  i  i  )  ^  ^  easily. 

Next  for  (ii)p,  and  ( iv) ^  ^  we  may  estimate  the  equality 


-10- 


ft  _l  r  _in-1,.0  ,  _in-1  ,  P  d”~%i.^)  . 


fo  r  f  0  )  .  D 


(L  )  •  D  p  dxdT 

X.  X, 


i  2y+y' 


ft  ,, 

/g  T  B  ( T)dT 


which  Implies  after  integrating  by  parts 


.  .  m  1 2  —  .  , t»t 

.  H  f  D  D  D  _m  i  _m 

t  J  +  — = —  Du  •Dp  dxj 

^  2^+ir'  ’‘i  't-0 


PP. 


2  ^m- 1 ^  ^m- 1 

>0  '  J  ^  ^  ^  ®x  *  °  Px 

2p+w'  2p+u'  i  i 


ft  I  r  P  i_m  .2  .  ^"^2 

L  T  /  — ^ —  Id  p  + - , 


_P _  Ir,"-!  i  |2.  . 

—  |D  I  dxdT 


2p+u' 


ft  .  £-1,.  ,0,t.  f  Id^pI  P  _m  i  _m  ,  . 

Jq  £t  (1-6  )  j  — ~ ^ —  D  u  •  D  p  dx  + 

2u+u'  ’'i 


+  C  ( t)t  +  B  ( t) t  dT  • 


These  equal ites  imply  (ii)„  and  (iv). 

"*  l,m 

Proof  of  Lemma  3.2.  For  any  positve  constant  e  >  0,  consider  the  form 


V  T  2  3 

(3.3)  y  (i)^  +  e  y  (ii)  +  e  (iii),  ,  +  e  (iv),  ,  + 

,  „  X  •  m  1,1  1,2 

)c=0  m=1 


By  taking  e  suitably  small  in  (3.3),  we  can  easily  prove  (3.2)  by  (3.3) 


Next  let  us  estimate  the  nonlinear  terms 


Then,  for  any 


Lemma  3.3.  Suppose  (p,u,9)  e  X(0,T>  E)  for  some  E  < 
positive  constant  e,  there  exists  a  positive  constant  C(e)  which  is 
independent  of  t  such  that 

2 

In  I  (1  +  T*')nA^(T)|  +  |B*'(t))  +  1C*'(T)I)  +  (1  +  T^)A^(T)dT 
0  »  « 

£=0 

<;  e  N^(0,t)  +  EC(e)N^(0,t)  . 

Before  proving  Lemma  3.3,  we  note  that  Lemmas  3.2  and  3.3  easily  imply 
the  desired  a  priori  estimates.  In  fact,  we  may  first  choose  e  so  small  and 
next  choose  E  so  small  that  we  have 

N^(0,t)  <  Cllpg,Ug,6Qll3  for  E  <  • 

Proof  of  Lemma  3.3.  Because  there  are  many  terms  to  estimate,  we  oick  up  some 
examples.  The  remaining  terms  will  be  estimated  in  the  same  way.  First  let 
us  pick  up 

(3.4)  J  f°p  dxdT 


in  Ip  A^(T)dT.  By  using  Nirenberg's  inequality  (2.4)  efficiently  we  estimate 
(3.4)  as  follows; 

I  1^  f  dxdil  =  I/q  [  - 

j  j 

=  '/J  / 

) 
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<  j  /q(supI pi ) api  iDundT 

X 

s  2 

<  C  /q  HpJ^  IlDull  ID^pl^dx 

1  i  1 

<  CE(aup«pl^  ]®(8up(  1+t)  I1DuI^)^[bup(  1+t)^IID^pI^)® 

T  T  T 

_  5 

X  /q  (1+t)  ^  dx 

<  CEN^(0,t)  . 


Next  let  US  consider  the  quantity 


(3.5) 


/„  /  D^(pu^)  .  D^p  dxdx 

0  X. 


rp^orlnq  In  t^A^(r)dx  or  T^B^(T)dT 


We  estimate  (3.5)  as  follows; 


<3‘5)  -  /n  /  T^(p  +  pu^ 


I-I..K-M-)  ’P  dxdx 

0  ^  ^  X  XXX  X  X 

j  jjcim  )i)<:m 


J-JxuDp  ‘Dp+Sxp  u  p  + 

•'0  •'  X,  x.x,  X  X  •  X  X  X 

j  ])t)lm  kHm 


2  j  2  3  j  3 

+  3xp  u  p  +xpDu  »Dp  + 

x.x,  X,  X  X,  x.x  ^x.  ^ 

3)cJlm)t)lm  3 


2  3  3  1  2  1 

+  X  (D  p  •  D  d)u  +3x0  u 


p  •  D  p)u  +  3x  p  u  p  + 

X,  X,  X  .x.x  X,  x.x 


^  2  I 

+  3x  p  u  p 


X,  X  .  X.X  X,  X . : 

K  3  m  K  j 


^  2  „3  j  ^3  ^  ^ 

+XpDu  ‘Dp  dxdx 
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<;  CE  Ig  <  CEN^(0,t)  , 


.  eN^(0,ti  --  fo  (/iD^l^d.fd  In^l^dxfaT 


N"(o.t>  ^  ""  !l  HD^un^  dT 


<  e.  N2(0,tt  (sup(1  +  t^«D^«^](3'iP(1^t^t(0^U«^^) 


X  t^t  +  T^'^dT 


<  e  N^(n,t)  +  Ce  ^EE^N^CO.t)  , 


,  IT^I  and  |lf,l  are  estimated  in  the  same  way 


as  1 1  - 


itj  -  ifs  f »' » •  ”’"1 


3  ..  „  4. 


<  I'"  /  sviplpl  «D  P"  “D 

I  3 


r  N  (n,t)  +  ^ 


‘M^n,t)  Jl  T  (I^T) 
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